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Abstract. Based on dynamical constraints about the Hubble se- 
quence evolution, observational data and a number of "conspir- 
acies", we propose that the dark matter around spiral galaxies 
OO is in the form of cold gas, essentially in molecular form and 
i— i rotationally supported. 

The existence of a much larger amount of cold gas in the 
(f-) outer discs of spirals is in accordance with dynamical studies 
taking into account the bar phenomenon and the secular growth 
of bulges, leading to the general conclusion that spiral galaxies 
have to evolve along the Hubble sequence from Sd to Sa over a 
couple of Gyr's. If so, the varying M/L ratio along the Hubble 
sequence suggests that dark matter is transformed into stars, 
i.e. dark matter should be in a sufficiently diluted form of gas. 

This hypothesis sheds light on several puzzles: 1) the con- 
spiracy of the flat rotation curves, 2) the estimated short gas 
consumption time of spiral galaxies, leading to the "gas con- 
sumption problem", 3) the constant ratio of dark matter to HI 
mass in the outer spiral discs, 4) the larger amount of visible 
gas in interacting galaxies with respect to isolated ones, giving 
rise to gigantic starbursts, 5) the high frequency of absorption 
lines in front of quasars, together with the large sizes of the 
absorbing medium around galaxies. 

Several problems raised by a completely gaseous dark mat- 
ter in the outer spiral discs are discussed. In particular the sta- 
bility and self-consistency of cold gaseous discs is shown to be 
a less severe problem than commonly believed. 
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1. Introduction 

Numerous assumptions have been made on the nature of dark 
matter, from modest extrapolation of well known objects, such 
as brown dwarfs, to a prolific variety of much more speculative 
hypotheses, such as WIMPS or modified gravity. Since dark 



matter seems to be required in different amounts over very dif- 
ferent scales, it is physically reasonable to expect that different 
kinds of explanations might be necessary along the cosmic hi- 
erarchy (Rees 1987). In some way the amount of dark matter 
is a measure of our ignorance, or of our lack of confidence in 
existing theories 1 , about the phenomena acting at each scale. A 
priori we have no reason to expect that our ignorance should be 
restricted to a single type of matter or physical effect. 

The present discussion is directed principally towards the 
problem of dark matter at the level of individual disc galaxies, 
for which it is the clearest. Therefore, a review about the main 
steps of the evolution of the ideas about dark matter inside these 
objects is useful. 

1.1. Dynamics in a historical perspective 

As well known, since Zwicky's early work (1933) dynami- 
cal evidence existed for large amounts of undetected matter in 
clusters of galaxies. Only forty years later, studies of field spi- 
ral kinematics by HI 21 cm observations (Rogstad & Shostak 
1972; Roberts & Rots 1973) gave a new impetus to this prob- 
lem at the scale of galaxies: flat rotation curves suggested that 
the presence of hidden forms of matter is a property of isolated 
spirals as well. 

In the early seventies two kinds of results appeared in con- 
junction: whereas the first flat rotation curves were published, 
seminal papers by Miller & Prendergast (1968) and Hohl & 
Hockney (1969) presented TV -body simulations of disc galaxies 
in view of explaining the formation and maintenance of spiral 
structures in discs. At that time, the galactic astronomers were 
focused on the linear density wave theory (Lin & Shu 1 964; Gol- 
dreich & Lynden-Bell 1965) concerning tightly wound spirals. 
The central problem for people working on TV -body numeri- 
cal experiments was to suppress these strong bars which form 
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1 For example most of the mass in the solar system, although 
hardly directly observable, is not considered as "dark" because 
the gravity of its bodies is consistent with the theories of their 
mass distribution. 
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spontaneously in the simulations, in order to obtain the well 
defined grand design spiral structure, "so often observed" and 
"well explained by the density wave theory". In fact, in the sev- 
enties, bars were not much studied and practically not observed, 
with the noticeable exception of de Vaucouleurs who insisted 
many times on their importance in the spiral classification. A 
further motivation to ignore them was that bars could hardly 
be included in the density wave theory. Essentially, two solu- 
tions were often considered to prevent the formation of bars: 
hot discs or massive halos (Ostriker & Peebles 1973; Einasto 
et al. 1974). But in the seventies no data could indicate that the 
inner parts of the galactic discs were hot enough. Discs were 
considered as consisting of cold populations by extrapolation of 
the solar neighbourhood data. Now, more recent observations 
indicate that the velocity dispersions in inner parts of typical 
Sbc galaxies can be of the order of 100 km s -1 (e.g. Lewis & 
Freeman 1989; Bottema 1993), showing that discs can be hot 
there. 

In the last decade, theoretical arguments, as well as new nu- 
merical experiments and deductions from the observations of 
the luminosity profile of discs, have contributed to modify the 
point of view about the role played by halos on global stability. 
The stabilisation of discs against non-axisymmetric perturba- 
tions, if ensured by halos, should only depend on the halo mass 
located in the inner region of the disc, where a bar usually devel- 
ops; it should not depend on the halo mass in the outer regions, 
where the rotation curve remains flat. Kalnajs (1987) disturbed 
minds by claiming 1) that the "luminous" rotation curves, cal- 
culated from the observed exponential stellar disc luminosity 
profiles, would fit very well the observed rotation curves in the 
optical region of galaxies, and, 2) that, as a consequence, dark 
halos would be useless to stabilise against a bar. 

Today, the global stability in the discs against bar-like per- 
turbations is no longer the acute problem that it used to be in the 
seventies: at least two thirds of observed spirals are recognised 
to be barred or to present an oval structure in the inner regions 
(de Vaucouleurs 1963). The remaining third include edge-on 
and dusty galaxies in which a bar can hardly be detected. In 
many cases, IR photometry, less affected by dust, reveals barred 
structures invisible in the B band; the most famous case being 
the Milky Way bar (Blitz et al. 1993). Further on, extensive nu- 
merical simulations by Athanassoula & Sellwood (1986) have 
shown that the growth of a bar may be prevented if the central 
part of the disc is hot enough. Finally, recent numerical experi- 
ments indicate how bars can be destroyed (Pfenniger & Norman 
1990; Pfenniger 1991a; Friedli & Benz 1993). Such dynamical 
constraints have implications on dark matter as explained in 
more detail in Section 2. 

1.2. HI rotation curves, maximum disc 

HI 21 cm observations of rotation curves have to be joined to 
the previous considerations. Spiral HI rotation curves, which 
extend typically up to 2— 3 -R25, are well known to stay in gen- 
eral sufficiently flat over this range for requiring an important 
presence of mass in the outer parts of spirals. Some authors tried 



to solve the question of the contribution of various components 
(disc, bulge, halo) to the observed rotation curves produced 
by HI data (see e.g. Sancisi & van Albada 1987). Among the 
proposed solutions, the "maximum disc" solution consists in 
maximising the contribution of the luminous matter to the ob- 
served rotation curve. This yields to a conservative lower limit 
for the amount of dark matter and an estimate of the M/L ra- 
tio for the disc. The maximum disc solution received recently 
a strong support from an extensive work on kinematical data 
of ~ 500 galaxies; it reproduces well the fine wiggles of the 
light, confirming that the fraction of dark halo within the stel- 
lar optical parts of the galaxies must be negligible (Freeman 
1992). Within the error margins, a small amount of dark matter 
in the optical disc is consistent with the estimates of the solar 
neighbourhood dark matter (Oort 1932; Bienayme et al. 1987; 
Kuijken & Gilmore 1989; Bahcall et al. 1992) which represents 
between and 50% of the local mass, though both extreme 
values are controversial. 

1.3. Conspiracy 

The only really obvious dark matter problem subsisting in spi- 
rals concerns the outer regions beyond the optical disc. The 
fact that the disc plus the bulge on one hand, and the dark halo 
on the other hand, bring an essential and about equal contribu- 
tion to the flat rotation curve in two distinct regions, namely in 
the inner luminous region and in the outer 21 cm emitting HI 
region, has suggested a physical coupling between these com- 
ponents, also called disc-halo conspiracy (Bahcall & Casertano 
1985; van Albada & Sancisi 1986). No convincing explanation 
of it has yet been given. 

Recently new HI observations have restricted the range of the 
conspiracy: Casertano & van Gorkom (1991) published HI ro- 
tation curves characterised by a large decrease between 1 and 
3 R25 and have found a clear correlation between the peak cir- 
cular velocity, its central brightness and the slope of the rotation 
curve in the outer parts of the discs. This result leads these au- 
thors to suggest that the ratio of the dark to luminous matter 
might be the critical parameter controlling the Hubble sequence. 
It is crucial to determine if really late-type galaxies contain more 
dark matter than early-type ones, or if the importance of dark 
matter decreases with luminosity. These questions introduced 
by Tinsley (1981) have been recently discussed again (see for 
example Salucci et al. 1991 and references therein). 

1.4. Bosma Hi/dark matter ratio 

A curious coincidence (Bosma 1981), until now unexplained, 
has been recently recalled by Freeman (1992, and references 
therein) and better documented by Broeils (1992): the surface 
density ratio of dark matter and HI gas in a sample of galaxies 
remains constant outside the optical disc, around 10-30. If 
dark matter tends to follow HI with a constant ratio of about 
20, we have the very curious situation that dark matter is con- 
comitant to HI radially, yet dark matter does not prevent HI 
to flare out the plane just beyond the optical disc. The flaring 
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shows at least that dark matter cannot be much flatter than HI. 
On the other hand, attempts to constraint the flattening of dark 
matter by the thickness of the flaring HI layer are either incon- 
clusive (Kundic et al. 1993), or indicate a flattened distribution, 
possibly as flat as HI (Rob Oiling, private communication, and 
Sect. 4.1). Certainly the simplest assumption is that if dark 
matter is proportional to HI radially, then it should also be 
proportional vertically. 

1.5. Overview 

To solve the previously summarised problems and conspiracies, 
we propose in this and a companion paper (Pfenniger & Combes 
1994, Paper II) the following hypothesis: dark matter in spiral 
galaxies is made essentially of cold molecular gas, in a mostly 
rotation supported disc. In the numerous conferences dedicated 
to the dark matter problem, extremely little critical attention has 
been devoted to this possibility. 

In Section 2 we review several new constraints that ad- 
vances in the large scale dynamics of galaxies have furnished. 
In particular the general conclusion that can be drawn from 
numerical dynamical studies is that spirals are unsteady objects 
and they have to evolve secularly along the Hubble sequence 
in the sense Sd — > SO, either by smooth evolution driven by a 
bar, or more dramatically by galaxy-galaxy collisions. Since the 
ratio of dark matter to stellar mass decreases along the Hubble 
sequence, from Sd to Sa, we can infer that stars should arise 
from dark matter, i.e. dark matter should be in a sufficiently 
diluted form of hydrogen, excluding Jupiters, brown and white 
dwarfs, and black-holes. 

In Section 3 we review the different observational con- 
straints on the possible forms of hydrogen. A massive amount 
of hot or warm form of hydrogen can be ruled out, essentially 
because hot and warm gas already fills most of the interstellar 
volume at a much too low density to contribute to the mass in 
an appreciable amount. On the contrary, although cold gas fills 
a small part of the volume, its density is large enough to encom- 
pass most of the ISM mass. Moreover, this condensed structure 
explains in part the invisible character of the medium, its low 
cross-section for absorption studies, and its transparency to ex- 
ternal radiation. Even if the medium is essentially molecular, 
the CO molecule can hardly serve as a tracer, because of the 
low metallicity of this quasi-primordial gas. A fraction could 
also be in atomic form, but HI emission cannot be detected in a 
medium at 3 K. Absorption studies of the outer parts of galaxies 
could probe this medium. 

Based on the fact that cold gas is observed to be fractal over 
several decades of length and density, a more specific discussion 
of the possible state of cold gas is presented in Paper II. It turns 
out that both the problem of mass underestimate in HI discs and 
the problem of star non-formation in outer discs is closely linked 
to the fractal structure. The physical state of this gas must be 
high density and cold temperature. Since no significant heating 
sources in the outer discs presumably exist, we can assume that 
the gas is bathing in the cosmological background, and that its 
temperature is about 3 K. In these nearly isothermal conditions, 



clouds can fragment until they reach small clump units, where 
the cooling time becomes comparable to the free-fall time. The 
average typical density of these elementary cloudlets, called 
"clumpuscules", is 10 9 cm~ 3 , column density 10 24 cm~ 2 , size 
30 AU, and mass 10~ 3 Mq . These small units are the building 
blocks of a fractal structure, that ranges upwards over 4 to 6 
orders of magnitude in scales. They are gravitationally bound, 
and the corresponding thermal width along the line of sight for 
molecular hydrogen at T = 3 K is about 0. 1 km s~ 

In Section 4 several problems raised by a cold gas dark 
matter are discussed, such as the stability and consistency of 
self-gravitating cold gas discs, the flatness of rotation curves, 
the effect of intense large scale magnetic fields, and implications 
at scales larger than galaxies. 

Several of the hypotheses and conclusions in this and 
the companion papers have a speculative character. Only ob- 
servations, particularly absorption line studies, will allow to 
test whether dark matter in galaxies can be, like in the solar 
neighbourhood, mostly or perhaps entirely baryonic, and, more 
specifically, in a form of hydrogen sufficiently diluted to permit 
star formation. 



2. Constraints from Theory and Simulations 

Starting in the late 60's, global dynamical studies of galaxies, 
in particular TV -body simulations, have profoundly changed the 
representation we should have of galaxies. Nowadays people 
familiar with TV -body simulations no longer consider galaxies 
as rigid objects formed for once 10— 15Gyr ago in a short 
time and remaining frozen until now, as proposed in the clas- 
sical ELS scenario (Eggen et al. 1962). On the contrary, discs 
are viewed as fragile structures (e.g. Quinn 1987) since they 
respond strongly to typical perturbers even as light as the Mag- 
ellanic Clouds. Interactions intermittently shake galaxies and 
can trigger important large-scale transformations. Tidal inter- 
actions have been extensively studied (e.g. Toomre & Toomre 
1972; Barnes 1988; Hernquist 1990). On the observational side, 
evidences for a decreasing but non-negligible amount of inter- 
actions and mergers with distance to galaxy cluster centers 
are now much better documented (e.g. Whitmore 1993). Also 
many cosmological indications, such as the variation of quasar 
density with redshifts, or the Butcher-Oemler effect (Butcher 
& Oemler 1978), indicate that galaxies and their environment 
have changed significantly within a time-scale of 10 Gyr. 

2.1. Disc formation and instabilities 

First, the formation of discs by gravitational collapse in a tur- 
bulent gaseous medium is a process badly characterised by a 
unique time-scale (as assumed in the ELS scenario), but by a 
range of widely different time-scales. Since Larson's (1969) and 
other's simulations of gaseous gravitational collapses with non- 
negligible pressure, the initial free-fall time Tff = (Gpo)~ ] / 2 is 
known to be the shortest time-scale of the collapse characteris- 
ing only the central region. In an indefinitely extended medium, 
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the collapse rapidly becomes inhomogeneous, pressure gradi- 
ents develop, and the collapse proceeds at larger radii for ar- 
bitrary long times. Since dynamical processes are much faster 
in the inner parts than in the outer parts, the inner parts can al- 
ready form a quasi-steady disc, while the outer parts still behave 
chaotically. For several tens of rotation periods the boundary 
conditions for the inner parts are not the ones of an isolated 
object: the concept of "formation epoch" is misleading. The 
slow growth of discs is vividly illustrated in high resolution 
cosmological TV -body simulations with gas (e.g. Evrard 1992). 

Second, dissipative discs are characterised by a wide choice 
of gravitational instabilities with relatively fast growth rates 
(pa lOGyr" 1 ) (see e.g. Binney & Tremaine 1987, Chap. 6). 
Owing to the fact that angular momentum is less easily dissi- 
pated than energy, a cooling disc tends to minimise kinetic pres- 
sure while conserving circular motion. Inevitably at some point 
the Safronov-Toomre stability parameter Q (Safronov 1960; 
Toomre 1964) has to become critical, i.e. the disc becomes 
gravitationally unstable, highly responsive to perturbations. 

A surprising large number of spirals, particularly the gas- 
rich Sc's, look chaotic, strongly suggesting a non-equilibrium 
state of these galaxies. Several TV -body experiments (Hohl 
1971; Sellwood & Carlberg 1984) display also chaotic spiral 
structures when an inelastic gas component is included. The ba- 
sic mechanism which appears to be able to regenerate the spiral 
structure in galaxies beyond a few revolutions is the so-called 
"swing amplification" (Toomre 1981). However such a process 
heats the discs rapidly to such an extent that it needs a constant 
supply of fresh cold gas to be maintained. Gas infall could be 
the solution of the problem of maintaining the chaotic spiral 
structure, as suggested by Toomre (1990). The privilege that 
only Sc's would receive gas seems strange unless earlier type 
galaxies would be former Sc's having evolved toward other 
types. 



2.2. The bar-bulge connection 

Typically, the evolution of a critically stable TV -body disc is 
that it easily forms a bar and spiral arms, allowing then a fast 
global transfer of angular momentum. Therefore the classical 
assumption that galactic discs remain almost unchanged over 
more than lOGyr is also easily invalidated. 

TV -body simulations have shown early (Hohl 1971) that a 
bar reorganises a stellar disc into precisely the observed charac- 
teristic exponential profile in the outer region (Freeman 1970) 
with a steeper profile in the bar-bulge region. At that time 
Hohl's result was not taken seriously into account since bars 
were viewed as anomalies. But now they appear instead to be the 
rule (a lucid explanation why bars are more natural than perfect 
axisymmetric discs in the region with a linearly raising rotation 
curve was given by von Weizs acker already in 1 95 1 , but forgot- 
ten during the following decades). More accurate simulations 
have confirmed Hohl's finding (Combes et. al 1990; Pfenniger 
& Friedli 1991). Therefore the dynamical effects of bars cannot 
be ignored. 



Another finding of TV -body work (Combes & Sanders 
1981), also not taken seriously first (Sellwood & Wilkinson 
1993), was that self-consistent bars excite important instabili- 
ties transverse to the galactic plane, tending to thicken the disc 
in z and generating a peanut-shaped bulge aspect of the bar. 
Bars are sufficiently non-axisymmetric structures for widening 
significantly all the orbital resonances, increasing chaotic mo- 
tions. In addition to the corotation, Lindblad and other radial 
resonances, bars have important vertical resonances able to am- 
plify star oscillations about the galactic plane (Pfenniger 1984). 
The variety of numerically generated peanut-shaped bars is re- 
stricted, and their density profiles quantitatively resemble the 
observed light profiles of peanut-shaped galactic bulges such as 
NGC 128 (Friedli & Pfenniger 1990). The connection between 
bars and peanut-shaped bulges was also confirmed by indepen- 
dent, more accurate, and more numerous calculations (Combes 
etal. 1990; Rahaetal. 1991). 

Recent studies of secular effects in discs have shown that 
the accumulation of gas, or any kind of matter, within the bar 
region (more precisely within its inner Lindblad resonance) 
is able first, to produce a secondary smaller bar (Friedli & 
Martinet 1993), and second, to later dissolve the bar into a hot 
spheroid half as large as the bar (pa the disc scale-length) akin 
to a small bulge (Hasan & Norman 1990; Pfenniger & Norman 
1990; Friedli & Pfenniger 1990). The required quantities of 
matter are only of the order of 2% of the galaxy stellar mass, 
illustrating the fragility of discs. In usual spirals such a mass 
transfer provoked by the gravity torques on the gas can be 
expected over a span of a few Gyr. So stellar bars appear to 
be fragile against dissipative perturbations and to be able to 
dissolve into small bulges especially fast in gas rich galaxies. A 
bar can eventually reform later if the disc cooling and sufficient 
gas accumulation bring the disc again to instability. After each 
bar dissolution the stellar mass of the bulge increases, since 
bulge building is an irreversible process. 

An alternative possibility to form bulges is accretion of 
satellites. However, the formation of small bulges such as Sc's 
bulges by accretion is a very unlikely process, because it re- 
quires to dissipate the energy of accreting satellites with an 
atypical small amount of angular momentum without heating 
the rest of the disc. In other words all the satellite kinetic energy 
must be deposited near the centre at the first passage, because 
inevitably subsequent passages cross the disc and heat it. In 
fact, TV -body simulations show that the repeated merging of 
dwarf satellite galaxies 1) destroys a bar (Pfenniger 1991a) and 
2) thickens a disc into possibly big bulges: 10% of such a kind 
of mass accretion can produce large bulges like in the Sombrero 
galaxy (Pfenniger 1992, 1993b). 

2.3. Energetics and disc sizes 

Simple global energetics considerations (Pfenniger 1991b, 
1992) show that the power emitted by stars in galaxies is com- 
parable to the gravitational power L gm that a galaxy can ex- 
change by large scale dynamics (L gm pa V s /G is the ratio of 
the gravitational energy GM 2 / R to the dynamical time R/V of 
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a virialised system RV 2 = GM, where V is the virial velocity, 
M the mass, R the radius, and G the gravitational constant). 
Thus, unless disc galaxies would be completely transparent, 
which are not since a substantial fraction of the light is ther- 
malised and reemitted in the (far) infrared, on the long run stellar 
light represents a power sufficiently large to counteract gravita- 
tion, leading to the dimensional relation L ho \ ~ V s /G linking 
the bolometric luminosity L\, \ to the gravitational power L gla . 
Comparing the actual luminosities and inserting the rotation ve- 
locities in this relation gives a striking match to the Tully-Fisher 
relation, in view of the invoked order of magnitude arguments. 
Note that the Lboi ~ V 5 /G relation applies only to a com- 
pletely self-regulated star forming disc. The lower exponent 
usually found in the IR Tully-Fisher relation can result from the 
contamination of young disc light by bulge and old disc light. 
To first order an old stellar population has a constant L/M ratio, 
therefore it should follow closely the virial theorem: L ~ V 2 R. 
The observed correlation for the ellipticals, called the "funda- 
mental plane", is indeed very close to this relation (Djorgovski 
1994). 

The power coincidence Lboi ~ L gri in actively star forming 
discs can only be understood if star formation is effectively self- 
regulated by a fast feed-back with global dynamics (Quirk 1972; 
Kennicutt 1989, 1990). It means that stellar activity 2 (internal 
physics), rather than initial conditions, determines galaxy sizes. 
Exactly the same is thought to be correct for stars, i.e. the 
relevant factor determining the size and mass of stars is not the 
initial fluctuation spectrum of the diffuse ISM, but the internal 
physics (mainly the atomic properties of hydrogen). 

If star formation in discs is self-regulated with the disc 
global dynamics, then it implies again that the aspect of disc 
galaxies must change over a time-scale comparable with the 
gas consumption time-scale, often found much shorter than the 
Hubble time with the conventionally determined amount of gas. 

2.4. Evolution along the Hubble sequence 

From the above considerations a more definite picture of the 
secular evolution of galaxies emerges. Disc galaxies can form 
from large scale and slowly rotating primordial gas condensa- 
tions as poorly structured, bulgeless and mostly gaseous discs, 
such as Im, Sm galaxies. As gas cools rapidly by radiation, 
proto-galaxies tend to flatten and to become more axisymmet- 
ric. Stellar activity is the principal energetic factor controlling 
the galaxy sizes. 

The virial theorem tells us that the energy loss increases the 
rotation speed simultaneously to the binding energy: slowly 
rotating Sd's are energetically less evolved than fast rotat- 
ing Sa's. The square of the typical rotation speed ratios 
((lOOkms-'/SSOkms" 1 ) 2 » 1/12) indicates the ratios of me- 
chanical energy dissipation. The increasing number of rotations 



Here "stellar activity", not to confuse with "chromospheric 
activity", is meant to include all the phenomena linked with star 
formation and evolution, and with all the mass, momentum and 
energy exchanges occurring between the stars and the ISM. 



progressively symmetrises the outer disc, and winds up the spi- 
ral arms. Effectively Sd's are typically asymmetric, Sc's have 
patchy and well open spiral arms, while Sa's-SO's are almost 
symmetric and have tightly wounded spiral arms. 

General cooling brings discs towards gravitational instabil- 
ity, leading to the formation of a bar, that secularly forms a 
bulge about half as large as the bar, so a small bulge. Beyond a 
critical accumulation of mass near the centre, bars dissolve into 
a spheroidal shape. Big bulges larger than the disc scale-length 
can result from the merging of a few dwarf satellite galaxies. 
So, secularly the mass of the bulge increases. Bulge growth is 
indeed an irreversible process; the Hubble sequence, from Sd to 
SO, is well known to be also a sequence of increasing bulge/disc 
ratio. 

The other irreversible process indicating the sense of evo- 
lution along the Hubble sequence is star formation (Kennicutt 
1990, 1992) since star formation continuously consumes gas 
and blocks most of it in long lived low mass stars. The conse- 
quence of nucleosynthesis, of stellar winds, and of supernovae 
is to increase the amount of dust and heavy elements. 

Therefore disc galaxies may only evolve along the Hubble 
sequence from Im-Sm-Sd to Sa-SO. The typical time-scale for 
a major change in a galaxy, such as the transformation from one 
Hubble type to the next, can be as short as 1 Gyr for dynamical 
factors and also of the order of 1 Gyr for star formation factors. 
But this time-scale is highly variable, since fortuitous merger 
events of arbitrary amplitude can boost not only dynamical 
evolution, but can trigger bursts of star formation. 

2.5. M/L ratio and transformation of dark matter into stars 

All these considerations about galaxy evolution have a direct 
implication for the nature of dark matter. If really galaxies do 
evolve along the Hubble sequence in the sense Sd to Sa, with 
a typical time-scale shorter than 10 Gyr, we have to understand 
why the ratio dark mass to stellar mass decreases systematically 
from Im-Sm-Sd to Sa-SO by a factor 10 2 while the ratio dark 
mass to HI mass remains constant, around 10-30 (Broeils 
1992). 

Also the approximate conspiracy connected with rotation 
curves has to be contemplated in time. Detailed studies of the 
age-metallicity and age-velocity relations in the solar neigh- 
bourhood have concluded towards a nearly constant star forma- 
tion rate (SFR) in the local disc during its evolution (Carlberg et 
al. 1985; Meusinger et al. 1991). A nearly steady state in SFR 
is also found in extragalactic studies (Kennicutt 1983, 1992; 
Donas et al. 1987). However Larson et al. (1980) had noticed 
that the average time-scale for consuming all the disc gas by 
star formation is of the order of 4 Gyr, in contradiction with 
the absence of evolution of the SFR; this problem is known as 
the "gas consumption problem" (cf. also Tinsley 1980). In fact, 
the constant SFR concerns mainly the middle spirals (from Sb 
to Sc) but not the later types and irregulars, affected by bursts 
(Scalo 1988), or the earlier types, that appear to have a declin- 
ing SFR accompanying their gas depletion. If the SFR is indeed 
steady in a large isolated spiral, this could be obtained by the 
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gradual extension of the diameter of its optical disc into re- 
gions with fresh gas, which are also dominated by dark matter. 
The mysterious conspiracy of flat rotation curves follows then 
naturally. 

Often fresh gas infall is invoked to explain the persistence 
of star formation (e.g. Larson et al. 1980). But discs are fragile; 
the amount of mass falling in from arbitrary directions and with 
arbitrary speeds is restricted for galaxies such as the Milky Way 
to a few percent of mass in order to conserve the flatness of discs 
(Quinn 1987; Toth & Ostriker 1992), a much too small amount 
for compensating the average rate of star formation. 

Preliminary calculations of the gas consumption time, 
quoted by Kennicutt (1992), indicate that the gas consump- 
tion problem can be less severe than previously estimated when 
the recycling of the gas is not assumed to be instantaneous, 
and a higher star mass loss is adopted. Delayed gas recycling 
and large stellar winds extend the star formation time-scale by 
a factor two to four, still shorter than the Hubble time. So the 
question remains to understand how the fraction of dark matter 
can be less in early type galaxies than in late type galaxies, 
while presumably early type galaxies were late type galaxies 
5- lOGyr ago. 

These evolutionary considerations and apparent problems 
suggest that galaxy evolution would seem much simpler if dark 
matter were in a form able to generate later stars, i.e. fresh 
diluted hydrogen and helium gas. This hypothesis excludes not 
only exotic particles, but also dense objects such as brown or 
white dwarfs, Jupiters, or black-holes. 

Since along the Hubble sequence the estimated amount of 
HI is systematically of the order of 1/10 of the dark mass, 
while the ratio of star mass to gas mass varies by two orders of 
magnitudes, it is clear that if a systematic error in the gas mass 
determination would be as large as a factor 10-30, essentially 
no dark matter would be required in discs (Casertano & van 
Albada 1990). 

Therefore, we are led to the idea that more gas than previ- 
ously thought, sufficiently diluted in order to be able to form 
stars later on, should be available in the outer discs. This gas 
would be mostly supported by rotation, shrinking in radius only 
slowly due to the difficulty to dissipate angular momentum. The 
energy dissipation has to be small since outer HI discs seem to 
persist for time-scales much longer than the HI cooling time 
without apparent corresponding energy input. In paper II a so- 
lution is proposed in which a very clumpy (fractal) structure of 
cold gas (T « 3 K) with quasi-isothermal conditions can con- 
siderably decrease the energy dissipation rate. Star formation, 
destroying the isothermal conditions in the optical disc would 
increase the dissipation rate at the transition from the optical to 
the HI disc. Gaseous dark matter would be consumed by star 
formation. The most likely possibility is that both gas would 
accrete slowly inwards, because dissipation cannot strictly van- 
ishes, and the optical disc would spread out with time, because 
the optical disc, made of dissipationless stars, can hardly shrink 
in radius. The exact amount of both inward gas motion and 
outward spreading of the stellar "fire" cannot be presently well 
estimated theoretically owing to our rudimentary knowledge of 



the involved physics of the ISM, star formation, and galactic 
dynamics. 

3. Observational Constraints on the Dark Gas 

We review in this Section the constraints that observations im- 
pose on the possible forms of mainly hydrogen. 

3.1. Halo gas in the Milky Way 

3.1.1. Halo gas in the solar neighbourhood 

The halo gas in our own Galaxy has been surveyed through 
optical and ultraviolet absorption studies in front of halo stars 
(e.g. Edgar & Savage 1989; Danly et al. 1992), and in front of 
remote quasars or X-ray emitters (e.g. Morton & Blades 1986; 
Kinney et al. 1991). An important result from UV studies is that 
the low ionisation material is about a factor of 10 times more 
abundant than the high ionisation material (e.g. Danly 1991). 
From the comparison of 21 cm HI emission towards the same 
line of sight as the absorbed halo stars, the z-distributionof halo 
gas can be derived. Danly et al. found evidence in most lines 
of sight for gas beyond z = 1 kpc, with column densities of the 
order of 1 — 12 • 10 19 atoms cm -2 . If the gas layer is fitted by 
an exponential z -profile, the exponential scale-length is 5 10 pc. 
There is an excess of negative velocities for the high-z gas seen 
in absorption, which has supported speculations on galactic 
fountains as the origin of this gas (e.g. Shapiro & Field 1976; 
Shapiro 1991). No evidence of neutral gas above the lower halo 
do exist: no further absorption lines occur towards extragalactic 
sources with respect to high latitude halo stars (above 2-3 kpc). 

Diffuse ionised gas has been also studied through pulsar 
dispersion measures. Since millisecond pulsars have been re- 
cently discovered in globular clusters at high z (~ 4 kpc), the z- 
distribution of the ionised medium has been derived (Reynolds 
1989). A characteristic scale-height of 1.5 kpc has been found 
for free electrons, significantly greater than that from neutral 
gas studies (Dickey & Lockman 1990). Highly ionised gas, 
traced by absorption and emission lines in UV, provides also 
evidence for a large exponential scale-height, roughly 3 kpc 
(Savage & Massa 1987). The volume filling factor of this com- 
ponent is low, between 3 and 10%. There is some indication for 
a substantial increase of the highly ionised species (Si IV, C IV) 
above 1 kpc (Pettini & West 1982), suggesting an extragalactic 
ionisation source. 

Above 1 kpc, the gaseous medium is therefore dominated 
by the diffuse ionised gas, that could be in transition regions 
between the HI clouds and the hot coronal gas. The source of 
ionisation, either from young stars or supernovae remnants from 
the plane (McKee & Ostriker 1977), or from the extragalactic 
EUV background (York 1982; Heisler & Ostriker 1988) is still 
subject of debate. The observed amount of NV appears to re- 
quire the existence of collisionally ionised gas near 2 • 10 5 K 
(Savage et al. 1990). 

Molecules have also been detected at high latitudes, in emis- 
sion (Magnani et al. 1985), and absorption (de Vries & van 
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Dishoeck 1988; Welty et al. 1989), but never at z larger than 
250 pc. Martin et al. (1990) have detected the H 2 fluorescence 
in the Lyman band towards the diffuse interstellar medium. 
They detect molecular hydrogen at high latitude, even when 
CO remains undetected because photodissociated, and suggest 
the presence of molecular halos, resulting from the photodisso- 
ciation of dense clumps. 

We can deduce from the Milky Way emission and absorption 
studies that there is indeed evidence for cool gas at high lati- 
tude; however, the possibility of disc-halo interactions (galactic 
fountains, bubbles, chimneys, etc.) prevents us to derive col- 
umn densities for the primordial halo gas. Also, if the dark gas 
around galaxies is a rotationally supported component, more 
in the shape of a thick disc than in the shape of a non-rotating 
spherical component, then we are at a unfavourable place in 
the Milky Way to detect such a component; external galaxies 
should give better constraints. 

3.1.2. Global disc and halo gas 

Only 7-ray measurements can inform us about the total number 
of nucleons integrated along the line of sight. They result from 
the interaction between cosmic ray protons and matter. The 
Galaxy is virtually transparent to 7-rays (Bloemen 1989), but 
the gas distribution is not known independently from the cosmic 
ray (CR) distribution. Therefore, no exact estimate of the total 
gas amount can be derived. In fact, the exponential scale-length 
of the 7-ray radial distribution is 15 kpc, much larger than that 
(~ 3 — 5 kpc) of the potential CR sources such as supernovae, 
pulsars, and massive stars, or also main sequence F to M stars 
which may contribute significantly to CR by sunspot activity 
(Casse & Goret 1978; Meyer 1985). This is interpreted in the 
frame of CR diffusion models (e.g. Bloemen 1989), but could 
also trace some hidden form of hydrogen. 

Indeed, first CR isotopic abundances studies show that the 
bulk of cosmic rays near the Sun comes from a distance of 
about 3 kpc from the galactic centre (Maeder & Meynet 1993). 
Second, CR particles are closely attached to the interstellar 
magnetic field lines, because of their small gyration radius, 
and the magnetic field intensity is a function of the volume 
density. The large radial extension of 7-rays point towards a 
larger volume density of gas. At large distances however, the 
gas volume density falls off at least as R~ 2 , implying also a 
fall-off of CR. Most of the outer dark gas could then not be 
traced by 7-rays since the amount of CR sources is low there. 

3.2. Gas around external galaxies 
3.2.1. Neutral gas 

In the 21 cm emission line neutral gas is detected in most galax- 
ies out to 1 —2Rh at a level of about 10 20 cm~ 2 , where Rr 
is the Holmberg radius (e.g. Bosma 1981; Sancisi 1988). Re- 
cent high-sensitivity observations have reached currently levels 
of 10 18 cm -2 , without increasing the HI sizes of galaxies. Van 
Gorkom et al. (1993) have undertaken a deep VLA survey of 
the galaxy NGC 3198 to test the HI extent at a resolution of 1 



arcmin (3 kpc), and discovered a considerably sharp decline of 
HI emission, suggesting the existence of an ionising front. The 
HI emission extends only exceptionally beyond 2R H (Hucht- 
meier & Richter 1982; Carignan & Freeman 1988). Many of 
these systems have an extended HI distribution due to a tidal 
origin (Weliachew et al. 1978; Simkin et al. 1987). An inter- 
esting feature in the outer parts of galaxies is the HI shoulder 
noticed by Sancisi (1988): after the optical disc radius, the radial 
distribution of HI surface density changes its slope (e.g. Wevers 
et al. 1986). The radial distribution varies then inversely pro- 
portional to radius, as the dark matter (Bosma 1981). The radial 
HI extent is also considerably reduced for galaxies in clusters, 
especially in the center of X-ray emitting clusters (Cayatte et 
al. 1990). This Hi-depletion is attributed to ram pressure strip- 
ping by the hot intracluster gas. 

3.2.2. Molecular gas 

The molecular gas radial extent is unknown in spiral galaxies. 
Only the CO emission yields indirect information on the H2 
content, and in general CO is observed only within the optical 
disc (e.g. Young & Scoville 1991, and references therein). In our 
own galaxy, weak CO emission has been detected at large radii 
(Wouterloot et al. 1990; Digel et al. 1993), somewhat outside the 
optical disc, but radial distances are hard to measure outside the 
solar circle. Recently, Allen & Lequeux (1993) observed also 
CO absorption in front of extragalactic continuum sources, at 
large radii. However, the CO/H 2 conversion ratio depends on the 
metallicity, and radial gradients of abundances are well known 
(Pagel & Edmunds 1981). Moreover, the surface density of the 
gas is rapidly decreasing with radius. If there exists a diffuse 
phase, the CO molecule looses its self-shielding well before 
molecular hydrogen, and could be entirely photodissociated, 
as in diffuse clouds in the solar vicinity (Blitz et al. 1990). 
Therefore, and essentially for metallicity reasons, the absence 
of detected CO emission does not preclude the existence of 
large amounts of H 2 . 

3.2.3. Diffuse ionised gas 

Diffuse ionised gas could be traced in external galaxies by 
Ha emission. For example Dettmar (1990) discovered a large 
extension (several kpc) of diffuse ionised gas perpendicular to 
the plane of NGC 891. He derives an electronic typical scale- 
height of 600 pc. However, this large z-extent is restricted to 
the inner half of the optical disc, and is highly related to star 
formation. Rand et al. (1990) founds an extended filamentary 
structure, with an even larger scale-height, between 2 and 3 kpc. 
The ionised gas layer is therefore thicker than what is obtained 
for the Milky Way by Reynolds (1989). This could be related 
to the higher star formation rate in NGC 891. The high-z Ha 
emission reveals a cut-off at R = 5-6 kpc. In the outer parts, Ha 
observations are not sensitive enough to detect column densities 
as high as 10 21 cm -2 , spread over a line of sight of 3 kpc, and the 
corresponding constraints on ionised gas around spiral galaxies 
are not significant (Braine & Combes 1993). 
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3.3. More molecular gas in interacting galaxies and 
starbursts 

Recently, Braine & Combes (1992, 1993) have shown that in- 
teracting galaxies, showing signs of tidal perturbations, have 
stronger and more centrally concentrated CO emission than un- 
perturbed systems. If the CO/H 2 conversion ratio is roughly 
constant, according to the non-varying molecular line ratios, 
then this implies a larger amount of molecular gas in interact- 
ing galaxies, by a factor 4-5 with respect to unperturbed ones. 
This larger amount is normalised to the blue luminosity, or to 
the physical size of the systems. The HI gas, on the contrary, is 
observed in about the same amount in interacting galaxies, with 
respect to unperturbed ones. Only for the stronger interactions 
and mergers is the HI content lower. 

Braine & Combes (1993) suggested a scenario in which 
large hydrogen reservoirs exist in the outer parts of galaxies, 
in rotational equilibrium. The strong gravity torques induced 
by an interacting companion, could drive progressively the gas 
towards the center (Combes 1991b). Gas infall could be more 
violent and sudden during the merging of equal-mass galaxies 
(Barnes & Hernquist 1991). Some evidence for such a sudden 
infall of "primordial" gas is found in the isotopic line ratios of 
molecular emission in mergers (Casoli et al. 1992). 

Fresh gas infall from the outer parts of galaxies helps to 
understand the existence of huge starbursts in mergers, and 
their exceptional infrared luminosities. In particular, a molec- 
ular mass of more than 10 11 Mq has recently been discovered 
in a distant infrared galaxy, of a total Lir = 10 14 Lq (Rowan- 
Robinson et al. 1991; Brown & van den Bout 1991; Solomon 
et al. 1992). The fact that the HI content begins to decrease in 
violent mergers, may suggest that they are running out of gas, 
or that the infall rate is insufficient. 

3.4. QSO absorption lines 

Maybe the clearer evidence for the existence of a large amount 
of gas around spiral galaxies is the widespread detection of ab- 
sorption lines in front of quasars. The large frequency of Ly-a 
absorptions (up to 100 systems in a single line of sight), remains 
unexplained. The derived mean cross-section for Ly-a is huge: 
assuming that they originate in galaxies puts the effective radius 
of spiral galaxies to 480 kpc, at a mean redshift of z = 2.5 (Sar- 
gent 1988). Also metal-enriched low-ionisation systems with 
TV (HI) larger than 10 18 cm -2 are found with large sizes, about 
3-4 R H (Bergeron & Stasinska 1986; Bergeron & Boisse 1991). 
Cases of 21 cm HI line absorbers of remote quasars by high- 
redshift systems have been reported (e.g. Briggs 1988). The 
derived column densities are of the order of 10 21 cm -2 . Two 
low redshift absorbers have been identified with nearby galax- 
ies (Haschick & Burke 1975; Carilli & van Gorkom 1987), and 
heavy elements are also present (Boksenberg & Sargent 1978; 
Rubin et al. 1982). However, Carilli et al. (1989) demonstrated 
that at least in one of these cases, the extended HI is due to a 
tidal tail, and proposed that most low redshift absorption line 
systems belong to interacting galaxies. 



Absorption systems with damped Ly-a lines are also char- 
acterised by their high column densities (larger than 10 20 cm~ 2 ), 
and are thought to be associated to galactic discs (Wolfe 1988). 
Their occurrence frequency along the line of sight of quasars is 
larger than normal spiral galaxies; this implies that spiral discs 
at that epoch (z ~ 2) had column densities of 2 • 10 20 cm -2 up 
to 3.5 R H , instead of 1.2 Rh measured in nearby spirals in HI 
emission, assuming the same comoving density of disc galaxies. 

As for molecular hydrogen, only one case, at z = 2.8, of 
H 2 absorption has been reported (Levshakov & Varshalovich 
1985; Foltz et al. 1988) with an estimated column density of 
~ 10 18 cm -2 . Observations are however difficult to perform 
and to interpret, given the multiplicity of lines in the Lyman 
and Werner bands, and the confusion with the Ly-a forest. 

Recently the Ly-a forest has been detected at low redshift 
with the Hubble Space Telescope towards 3C273 and PKS 
2155-304 (Morris et al. 1991; Bruhweiler et al. 1993). The 
large number of lines detected (9—16) was a surprise; the Ly-a 
absorbers being much more frequent than expected by extrap- 
olation of the high-redshift evolution. The number of lines was 
even compatible with no evolution since z = 2. According to 
their metallicity and clustering characteristics, the absorbers ap- 
pear to have the same basic properties as those at high redshift. 
The abundance of gas around galaxies is therefore maintained 
at low redshift. An alternative interpretation has been recently 
proposed in terms of ionised gas (Maloney 1992), extending a 
few hundreds of kpc around galaxies. 

3.5. Assumptions leading to neutral gas mass 

When the hydrogen mass is determined from the 21 cm line ob- 
servations, several assumptions are made, the violation of each 
one leading to an underestimate of the mass. Recent extensive 
reviews about HI observations have been written by Kulkarni 
& Heiles (1987) and Burton (1992). 

In galactic conditions the 21 cm line emission is almost 
always collisionally excited (for a density larger than about 
5 • 10 _3 cm -3 ). The actually measured signal is the difference 
between the brightness temperature and the background tem- 
perature, which is proportional to the difference between the 
gas kinetic and background radiation temperatures. Often the 
background temperature is neglected, but it is at least at 3 K. 
So any gas cold enough to have a temperature close to the 
background temperature is barely detectable, independent of 
the optical depth. Hence one should be careful before ruling out 
very cold conditions especially because the cold gas cooling 
time by black-body radiation is always very short with respect 
to the Hubble time. In fact, observations of CO and continuum 
radiation at 20 cm in the Galaxy and nearby galaxies suggest 
that if the main cold gas heating source is cosmic rays, then 
large amounts of unseen cold gas near 3 K should exist (Allen 
1993). 

A very important assumption usually made is that HI is 
optically thin. Extremely thin column densities below about 
10 18 cm~ 2 are difficult to detect, and for temperatures below 
100 K, column densities higher than about 10 21 cm -2 are opti- 
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cally thick. Therefore the range of measurable column densities 
in HI is practically less than 10 3 , and more typically of the order 
of 10 2 . Similar restrictions apply to other single line observa- 
tions (cf. Scalo 1990). Larger ranges of column densities can be 
investigated only by a series of lines. A priori the ISM possesses 
a range of densities between the intergalactic medium density 
(< 10" 3 Hcm- 3 ) and the stellar density (~ 10 24 Hcm" 3 ). In 
comparison, the highest measured densities below the ones of 
stars are of the order of 10 7 — 10 9 H cm -3 only (i.e. Mezger et 
al. 1988; Migenes et al. 1989; Wilking et al. 1989). These are 
condensations associated with star formation, i.e. warm cores. 
In some cases they correspond to accretion discs around stars. 

Observational surveys show that HI in the Galaxy can be 
both very transparent, since one observes well at large distances 
throughout the disc, and completely opaque, since the bright- 
ness toward the center and anti-center of the Galaxy are equal 
(as explained by e.g. Burton 1992, Chap. 2.4). Such seemingly 
contradictory conclusions are only possible if the HI distribu- 
tion is extremely clumpy. Despite the Doppler shift that makes 
HI more transparent to the 21 cm line, in several other directions 
at low Galactic latitudes the 21 cm optical depth hovers around 
t = 1 (Burton 1992). Therefore the assumption of general opti- 
cal thinness for the 21 cm line is not supported by observations, 
particularly because HI is known to be clumpy up to the high- 
est resolution. The main argument for optical thinness, that the 
21 cm line can be detected at large distances in the Galactic 
disc, is invalid if HI is far from being homogeneous, actually 
the general evidence. 

As for the molecular phase traced by the CO emission, it 
is well known that the CO line is generally optically thick; the 
CO/H 2 conversion factor is empirical, justified by the virial 
theorem (e.g. Combes 1991a). This ratio is unlikely to be valid 
when the metallicity is far from solar, or when the excitation is 
insufficient (Allen & Lequeux 1993). 



3.6. Formation of molecular hydrogen 

Molecular hydrogen is the lowest energy state of neutral hy- 
drogen. Naively one should therefore expect its frequent occur- 
rence in space. Dissociating UV radiation and slow chemical 
reaction rate in a low density medium around 1 cm -3 are usu- 
ally invoked to rule out large fractions of molecular hydrogen 
in the outer disc. 

Both assumptions can be vitiated in a hierarchical medium 
that both shields itself and contains dense clumps, and is remote 
from massive star formation sites such as in the outer galactic 
discs. Indeed the time-scale of H 2 formation by simple neutral- 
neutral reaction, 3 x 10 6 / «h Gyr (see e.g. Genzel 1992) is much 
longer than the Hubble time in a homogeneous medium with 
n H ~ 1 cm -3 , but becomes short with respect to galaxy ages in 
high density clumps with n H > 10 5 cm -3 . If ionising radiation 
or dust are present in some small amount, catalytic reactions can 
form H2 much faster than the neutral-neutral reaction. There- 
fore the above H 2 formation time-scale is a conservative upper 
bound. 



Such high density clumps are common in molecular clouds 
where molecules with heavy elements allow us to observe them, 
but molecules such as CO just trace mass and are unlikely the 
unique cause of these high densities. No CO/H 2 relation can 
obviously exist in a high density pure H+He medium in which 
no stars have ever enriched it with heavy elements. 

3.7. Fractal geometry of cold gas 

The fractal nature of the interstellar cold gas has been increas- 
ingly well documented (Falgarone et al. 1992). Nearby molec- 
ular clouds are observed self-similar in projection over a range 
of scales and densities of at least 10 4 , but perhaps up to 10 6 . 

According to measures of HI in VLBI and of variable emis- 
sion of QSO's, structures on a scale as small as 25 AU (Diamond 
et al. 1989), or 10 AU (Fiedler et al. 1987) have been detected. 
Church well et al. (1987) also have detected 22 compact ra- 
dio sources in Orion, whose sizes range from less than 27 AU 
to 230 AU. Cometary globules in the Helix planetary nebula 
(Meaburn et al. 1992) are examples of dense (> 6 • 10 5 H cm -3 ) 
and small (« 60 AU) clumps that seem to be evaporated by the 
central hot star. 

Remarkably the fractal dimension of cold gas isophotes at 
different wavelengths is almost constant, around 1.3—1.4, over 
a wide range of scales, and independent of the molecular or 
atomic state of the gas, or whether it is thought to be self- 
gravitating or not (Scalo 1990; Falgarone 1992). 

Modelling the physics of fractals is presently far from triv- 
ial, as one of the most basic traditional assumptions in physics, 
differentiability, is lost in practice, due to the huge dynami- 
cal range needed to describe the strong gradients. Unless the 
smallest scale of a hierarchical structure is accessible (to ob- 
servations, respectively to calculations), where differentiability 
might be a valid assumption, in principle the use of differen- 
tial equations is not justified if the resolution is insufficient. In 
particular the usual radiation transfer equation and the hydro- 
dynamical description need a differentiable medium. 

Often physical properties can greatly differ between smooth 
and fractal objects, because macroscopic properties may depend 
on the smallest scale at which the fractal behaviour stops. Al- 
though the general evidence is that HI is structured down to the 
smallest accessible scales, lacking better tools observers apply 
the equation of transfer and assume homogeneity at subresolu- 
tion scales to estimate the amount of interstellar gas, which in 
turn provides the amount of baryons in the outer discs of spirals. 

In a companion paper (Paper II), and in Pfenniger 1993a, 
we show by fractal models that large underestimates of the 
mass can result by factors more than 10 when assuming that 
the fractal dimension is 3 (diffuse medium) while in reality it 
would be less, around 1.5-2. 

3.8. Summary 

Only a small fraction of the mass can exist in a diffuse, partially 
ionised warm gas which would then correspond to the metal- 
lic lines of ionised systems observed in absorption in front of 
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quasars. This phase is to be expected, since the outer parts of 
galaxies are reached by UV and soft X-ray photons from quasars 
(Heisler & Ostriker 1988; Hasinger et al. 1991). Possibly the 
sharp HI emission fall-off in the outer parts observed by van 
Gorkom et al. (1993) is due to the ionisation of this diffuse 
phase. Since the diffuse phase fills most of the volume at a low 
density, it cannot contribute significantly to dark matter. 

On the other hand the existence of a cold, dense, self- 
shielded and very low volume-filling gas in fractal form is 
compatible with observations. This is the only form where a 
large amount of baryonic mass can be contained. A factor 10 
or more of hydrogen mass underestimate is enough to remove 
the need of exotic matter in disc galaxies (cf. Casertano & 
van Albada 1990). In addition to large optical depths, both the 
formation of H2 and a low temperature close to 3 K may con- 
siderably increase the amount of mass not emitting the 21 cm 
line. 



4. Discussion 

The subject of dark matter, even when restricted to spirals, cor- 
responds to such an extended literature that all its aspects cannot 
be discussed in a single paper. There has been good reviews on 
the dark matter problem with different emphases, by Trimble 
(1987), Carr (1990), Ashman (1992), and Tremaine (1992). We 
briefly mention below some points that seem relevant to us. 

4.1. Disc stability and self-consistency 

A question that arises immediately if disc dark matter is as- 
sumed to be made of cold gas is to explain the stability of such 
discs when the radial velocity dispersion is inferred from the 
observed HI or CO velocity dispersion (<r gas = 5 — 14kms _1 ). 
Linked to this is the problem of explaining how dark matter 
might be confined to a disc as thin as the flaring HI discs when 
considering the low velocity dispersion of the gas. 

In order to prevent gravitational radial instabilities the 
Safronov-Toomre stability parameter 



Q 



ttGS ' 



(1) 



(Safronov 1960; Toomre 1964; Goldreich & Lynden-Bell 1965) 
must satisfy Q > 1 . As usual k is the radial epicyclic frequency 
and £ is the disc surface density. This criterion is often evaluated 
with the Mestel's thin disc model £ ~ R~ 1 . In a such thin disc 
with flat rotation curve v c = constant, we have k = \/2v c /R 
and £ = v*/2irGR (cf. Binney & Tremaine 1987, Chap. 2.6). 
With v c = 200 km s -1 stability requires <r gas ps Qv c /t/S > 
7 1 km s~ 1 , independent of radius. 

In outer parts of disc galaxies <r gas is indeed observed to be 
almost constant (e.g. NGC 1058 by Dickey et al. 1990), which is 
consistent with a Q-regulated gaseous disc with a radial density 
profile proportional to dark matter (consistent with the Bosma 
Hi/dark matter ratio) . However with respect to 7 1 km s ~ 1 , <r gas 
is too small by a factor ps 10, which is usually taken as an 
argument for invoking an external thick dark matter halo. 



Actually observed HI discs appear as far from being per- 
fectly axisymmetric and smooth; they can present large asym- 
metries, warps, spiral arms and massive HI complexes, so the 
necessity of strict stability is not demonstrated (see e.g. M33 
and M101 as reviewed by Brinks 1990, or even the prototypical 
dark matter spiral NGC 3198, Begeman 1989). The irregular 
shapes of HI discs from large scale down to the beam reso- 
lution suggest on the contrary that the gas is unsteady, and a 
subcritical Q might be acceptable. Also when a disc is substan- 
tially non-axisymmetric, the applicability of Eq. (1) is less obvi- 
ous. Instead one might consider average quantities over several 
turns, or over 2ir in azimuth at a given time. Then the effective 
radial velocity dispersion is increased by the contribution of 
the radial non-circular motions 3 . To first order for a constant 
circular velocity v c , in the average disc the epicyclic velocity 
is ps kAR = \/2v c AR/R, where AR is the non-axisymmetric 
offset at the radius R. The radial component of the epicyclic 
velocity is smaller by a factor 1 /a/2. So the effective radial 
velocity dispersion is given by 



; eff 



vi + at, 



(2) 



For example for AR/R = l/\0, v c = 200kms _1 and <r gas = 
lOkms -1 , we get <r e ff ps 22kms _1 . 

Furthermore the HI discs are far from being thin. On the 
contrary they often flare out approximately linearly with radius 
beyond the optical disc, before being warped at larger distances 
(see e.g. Burton 1992, Merrifield 1992), giving thus the im- 
pression of an exponential flaring. The general effect of a finite 
thickness on the stability is to increase the effective Q by a 
factor 1 + kh z (Vandervoort 1970, Romeo 1992), where k is 
the wave-number and h z the disc scale-height, because mass 
being diluted in z the responsiveness of the disc is decreased. 
If we replace £ in Eq. (1) by 2p(R)h z , were p(R) is the mass 
density in the plane, and for flat rotation curves we have also 
p(R) = p(Ro)Rl/R 2 , we obtain 



h z = 



V2 ( a 



gas v c 
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So if the ratio within parentheses is constant, h z increases in- 
deed linearly with radius, which is precisely the behaviour ob- 
served by Merrifield (1992) in the Milky- Way for R < R < 
2Rq (see his Fig. [7]). Let us estimate Q for the Galaxy with data 
from Merrifield's paper: <Tg as = lOkms" 1 , v c = 200^8"', 
p(Ro) = O.lM pc~ 3 , R = 8kpc, and his Fig. [7] gives 
h z ps 0.045 R. Thus from Eq. (3) we extract Q ps 0.36. If 
a part of the observed flaring is contaminated by the warp, the 
true h z /R would be smaller and Q would be increased. If taking 
into account non-axisymmetric motions <7 ef f = 20 -30 km s -1 , 



3 It is easy to show that the continuity and momentum equa- 
tions have the same forms for time-averaged quantities, except 
that the pressure term is increased by the square of the time- 
average velocity variations. It follows that the same stability 
criterion can be applied for the average quantities. 
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we reach a domain of subcritical Q ~ 0.8 - 0.9. Finally, if we 
take the finite thickness into account, the lowest wave-number 
k being 2-k/R, we must multiply Q by a factor of at least 
1 + 2irh z /R 1.28, which yields a marginally stable disc, 
Q e ff 1. 

Let us check whether the above data is consistent with a 
completely self-gravitating disc as thin as the observed Milky- 
Way gaseous disc. For flat rotation curves we have dK z (z) /dz = 
4-irGp(z), where K z is the force along z. Integrating the Jeans 
equation 4 



P<T 2 Z (Z) 



pKz dz 



4ttG 



dK, 
dz 



Kz dz , 



(4) 



gives <r 2 (0) = K z (oo) 2 /8irGp(0). Approximating dK z /dz(Q) 
by K z (oo)/h z gives h z = a Z (Q) / '^2%G p(Q). Self-consistency 
with flat rotation curves requires p = p(Rq)(Rq / R) 2 , which 
yields to 



h z = 



t z (0) 



R 



V2vGp(Ro) Ro ' 



(5) 



Again h z is proportional to R if a z is constant. Assuming 
a z (0) = A/2cr gas (averaged over z, a z is certainly smaller 
than a z {0)) we get, using the previous numerical values, 
h z = 0.034.R, so quite close to the observed (and probably 
warp-contaminated) h z = 0.045.R taken from Merrifield's pa- 
per. 

Let us discuss the possible bias due to the effect of an 
outer disc flaring and a truncation at a finite radius on the 
usual mass estimate by the rotation curve. Indeed in a truncated 
Mestel disc the rotation speed increases inside and near the 
truncation radius, which obviously results in an overestimate 
of the inner mass. In several disc mass models with either 
constant height or linearly increasing height, we have compared, 
by numerical quadratures, the effective inner mass with the 
mass deduced from the virial relation GM = v 2 R. In order 
to conserve an almost flat rotation curve some mass models 
have been taken as p ~ arccos(i?/i? max )/i? 2 , where i? max is 
the truncation radius. The general result is that at R < R milx 
the inner mass is overestimated, with respect to the virial value, 
from k, 0% at R = up to < 30 - 40% at R = E max . An 
overestimate by 40% is certainly an upper bound for realistic 
discs with a mild truncation. The effect of flaring has been found 
negligible. Thus the disc truncation effect can only decrease the 
required dark matter by a small amount. 

Therefore taking into account the non-axisymmetric aspect 
of HI discs and their finite and flaring thickness alleviates the 
need of large <r gas . A <r gas of the order of ~ lOkrns -1 , ap- 
pears much less a problem in a realistic non-axisymmetric HI 
disc than in the perfect Mestel disc. Also linearly flaring and 
completely self-gravitating gas disc are possible with a velocity 
dispersion of the order of only 10 km s~ 1 . 



4 A reasonable assumption is that K z tends asymptotically 
towards a constant or a maximum at a z where p is negligible. 
Here the symbol oo is understood as the height at which K z 
reaches the maximum. 



In fact the problem of stability is not restricted to large scale. 
Cold gas in outer discs is strongly Jeans unstable already at small 
scales (< 100 pc), as clearly visible in the Milky Way by the 
profusion of details in the HI (e.g. Burton 1992), yet it survives 
there for Gyr without apparent energy input. In Paper II we ar- 
gue that the paradox of cold gas persistence without formation 
of dense objects such as Jupiters or stars can be understood by 
its fractal structure at a temperature close to 3 K. At the small- 
est scale of the fractal hierarchy the smallest clumps, called 
clumpuscules, are prevented to collapse by collisions because 
the collision time is much smaller than the collapse time. The 
smallest scale of the fractal is determined by the transition from 
the isothermal to the adiabatic regime. As a consequence of col- 
lisions the clumpuscules exchange energy with their neighbours 
and with the 3 K radiation. The near isothermality suppresses 
energy dissipation. The effect of the inclusion of the 3 K radia- 
tion for the energetics of outer discs does not seem to have been 
investigated. Indeed, in the outer discs the energy density of the 
3K radiation (u = 4aT 4 /c 0.26 eV cm -3 ) becomes com- 
parable to the turbulent kinetic energy of the cold gas beyond 
R « 20kpc. 



4.2. Flat rotation curves and global energy minimum 

Traditionally the flatness of the rotation curves is explained by 
the presence of a massive and hot collisionless isothermal dark 
halo at the kinetic temperature of the rotation speed. Now, if 
dark matter is in a cold gaseous form and self-gravitating, one 
should understand why the rotation curves tend to be flat. 

If most of the mass is cold, it has to be supported mostly 
by rotation and the condition of isothermality no longer leads a 
priori to constant rotation curves. However, the total energy E 
of a self-gravitating and rapidly rotating disc is easily seen to 
be a functional E[$] of the potential $ alone 



= / / 2%RdRdz {\pv 2 c +\p§) , 



where 



P 



4ttG 



and 



2 

c = ~dR 



(6) 



Minimising E[$] by variational calculus restricts the class of 
potential functions $ precisely to the ones having a flat rotation 
curve (Pfenniger 1989): 



SE[$] = requires = 

oR 



(7) 



This result is obtained without constraining the boundary con- 
ditions, just what is needed for slowly growing discs with time- 
dependent and chaotic outskirts. This result is robust in the sense 
that small perturbations to the total energy integral in Eq. (6) 
leads to small perturbations of the result in Eq. (7). Also pertur- 
bations produced by mass distributions outside the considered 
volume do not change the result. In any case, flat rotation curves 
are incompatible with boundary conditions of isolated objects, 
since mass diverges proportionally to the radius. 
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4.3. Magnetic fields 

Magnetic fields are sometimes invoked for solving the dark 
matter problem in discs (e.g. Nelson 1988; Battaner et al. 1992), 
with the implicit hope that by complicating the physics new 
alternatives can emerge. 

Yet an elementary argument 5 based on the virial theorem 
shows that any magnetic field increases the dark matter prob- 
lem. In the virial theorem (cf. Chandrasekhar 1968, p. 581) the 
different interacting energies of a self-gravitating systems in 
uniform expansion or contraction (and in particular in equilib- 
rium), sum up to zero. Except the gravitational energy, in an 
isolated system all the conceivable energies contribute to the 
sum positively, including the magnetic energy. 

Admittedly, if the system is non-isolated the magnetic (or 
gas) pressure at the boundary contributes negatively to the virial 
equation. But this negative contribution can outweight its inter- 
nal positive contribution only if the outer pressure is larger than 
the volume-averaged internal pressure. If they are equal, their 
effect cancel. 

Therefore increasing the magnetic field intensity inside a 
galaxy implies, that at the admitted known rotation speed and 
size of the galaxy more mass is required to compensate the 
additional magnetic pressure. The observed rotation velocity 
implies a lower bound on dark mass. Additional non-circular 
motions, gas pressure or magnetic fields increase the positive 
terms in the virial equation and therefore increase the required 
mass. 

Finally, if a continuous conducting gas can eventually be 
supported by magnetic fields at the galactic scale, this assump- 
tion is harder to justify in a highly hierarchical and turbulent gas 
in which densities vary by more than five orders of magnitude 6 . 
In any case young stellar populations born from the gas cer- 
tainly cannot be supported by a general magnetic field at the 
galactic scale. Therefore if gas would be supported mainly by 
magnetic fields and would rotate much faster than gravitation 
would require alone, then the orbital angular momentum of 
young stars would exceed the one of circular orbits, and, con- 
trary to observations, the velocity dispersion of young stellar 
populations would be much larger than the velocity dispersion 
of the gas. 

Other arguments against a significant role of magnetic fields 
have been presented by Persic & Salucci (1993): the kinematics 
of galaxy pairs, satellite galaxies, and small groups all require 
dark halos, not magnetic fields. 

4.4. Polar rings and warps 

As for galaxy companions, polar rings are indicators of the local 
gravitational potential, so they constraint the total inner mass, 
but they weakly constraint its distribution. 

5 During the refereeing phase of this paper very similar argu- 
ments have been published in Nat. 356 (1993) by J.R. Jokipii 

6 E.H. Levy (p. 19), and by P. Cuddeford & J. Binney (p. 20). 
6 In comparison the density difference between common gas 
and solid objects on Earth spans only 3-4 decades. 



Over the years polar ring studies have less and less con- 
strained the dark halo to be spherical. Recent works (Arnaboldi 
1992) indicate the possibility of relatively flat (c/a 0.6) ha- 
los. We remark that constraining the shape of the mass distribu- 
tion by the shape and kinematics of the orbits of test particles is 
an intrinsically difficult task. Major modifications of the galaxy 
mass distribution far from the polar ring change the orbit shapes 
very little, while, on the contrary, the self-gravity of a polar ring 
as light as a few percents of the total mass can not only deter- 
mine its stability, but also modify the orbit shape. This is a 
crucial point when the polar ring is considered as a test body in 
the potential of a dark halo (Sparke 1986). 

Concerning warps, it is interesting to recall an often forgot- 
ten paper by Petrou (1980). She showed that precisely a dark 
halo that becomes more flattened with radius may be a solu- 
tion to the problem of explaining why observed warps have a 
straight line of nodes. A priori the warp line of nodes should 
rapidly wind up in a differentially rotating disc. In the context 
of spherical, or slightly triaxial massive halos no convincing 
solution to the line of nodes problem has been proposed yet. 

A related argument for a flat dark halo in the outer discs has 
been put forward by Sancisi & van Albada (1987): If the outer 
gas in spirals may have been accreted recently, as suggested by 
its asymmetric and irregular shape, the fact that it has already 
almost found the optical disc plane indicates the existence of a 
massive flat disc in the outer parts, more or less coplanar with 
the optical disc. 

4.5. Dark matter at larger scales 

Kinematical studies concerning the galaxy rotation curves in 
clusters are still based on too small samples to infer general 
conclusions on the presence, or absence, of an important mass 
in the outer part of galaxies in such environments. More and 
better resolved optical as well as 21 cm data are necessary. 

A clear HI deficiency is observed for the galaxies in the 
inner parts of the clusters with respect to the more external 
regions, generally attributed to a stripping process. Whitmore 
et al. (1988) claim that in these galaxies optical rotation curves 
drop somewhat in the regions beyond the HI limit, while Amram 
et al. (1993) contest the reality of the declining rotation curves. 
It is too early to say if there is a correlation between the HI 
deficiency and a decrease of the rotation curve. 

At the scale of clusters, more dark matter appears to be 
present, with a density parameter in the range = 0.1 —0.3 
(e.g. Peebles 1 984). The possibility for this matter to be baryonic 
too is not excluded. Although this remains to be examined in 
details, some recent discoveries about hot gas in clusters makes 
it attractive. Mushotzky (1993) presented X-ray observations 
of galaxy clusters in which the hot gas of large clusters is able 
to represent a mass comparable to the dark mass. 

The distribution of dark matter inferred from the observa- 
tions of gravitational images in the shape of arcs is clearly more 
peaked in the center of the cluster than the visible matter (Ham- 
mer 1 99 1 ). This is also confirmed by the interpretation of the hot 
gas X-ray emission (Gerbal et al. 1992). In rich clusters, where 



D. Pfenniger et al.: Is dark matter in spiral galaxies cold gas? I. 



13 



most of the galaxies are of early type, so are not surrounded by 
a lot of cold dark gas, most of the gas must be in a common 
envelope, with a characteristic velocity dispersion which can 
be, at some scales, lower than the intra-cluster dispersion. 

Nucleosynthesis constraints permit the possibility of bary- 
onic dark matter, since the baryonic mass contribution to the 
necessary mass closing the Universe, O b in the frame of the stan- 
dard Big Bang nucleosynthesis could be as high as Ob = 0.1 
if the Hubble constant is within its lowest possible values 
(e.g. Boesgaard & Steigman 1985; Smith et al. 1993). However 
other uncertainties already in the frame of this model could raise 
Ob to 0.2 (Schramm 1991). In the frame of inhomogeneous cos- 
mologies O b can reach 1 (e.g. Malaney & Fowler 1988). The 
visible baryons already account for O b « 0.01, and we propose 
that around galaxies, the dark baryons increase this contribu- 
tion to O b « 0. 1 . In any case, the nucleosynthesis lower limit 
already requires a lot of dark baryons. 



5. Conclusions 

We have recalled that the arguments for the presence of dark 
matter in spirals have become weaker for the optical parts of 
galaxies, while dark matter is really needed in outer gaseous 
discs. Other pieces of information now suggest that disc galax- 
ies evolve faster than thought 30 years ago. In a couple of Gyr, 
a bar can form and dissolve, rearranging the disc into an ex- 
ponential profile, and heating the central disc vertically into a 
bulge. Already internal factors such as star formation and large 
scale dynamical instabilities speed up evolution, but external 
perturbations, such as satellite interactions and mergers, trig- 
ger also large scale instabilities and star formation. The rapid 
transformation of gas into stars, when contemplated with the 
dark to stellar mass conspiracy in the rotation curves, particu- 
larly points to look for the possible identification between dark 
matter and a form of matter able to form stars, i.e. sufficiently 
diluted gas. 

Better observations of the ISM show that the cold gas is 
fractal and essentially clumpy down to very small scales, of 
the order of a few tens of AU. A wide range of gas densities 
(> 10 9 H cm -3 ) that continuity between dense clumps and stars 
leads to infer the existence, are difficult to investigate by today's 
techniques. As shown in Paper II, HI mass determination could 
underestimate the gas mass by a factor 10 or more owing to the 
very inhomogeneous nature of cold gas leading to high optical 
depths in a small fraction of the sky, rapid H 2 formation, and 
thermal equilibrium with the 3 K radiation. The dark matter 
to HI constant ratio would seem then natural. If such a large 
error in the gas mass determination can be confirmed by new 
observations, then the problem of dark matter in discs would 
be solved. 
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